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We describe the POWHEG-BQX package, a general computer code framework for implementing 
NLO calculations in Shower Monte Carlo programs according to the POWHEG method. 



1 The POWHEG method 

Next-to-leading order (NLO) perturbativc QCD computations as well as Shower Monte Carlo 
(SMC) programs are fundamental tools for the present-days particle physics phenomenology. In 
particular, SMC programs incorporate the description of a generic high-energy hadronic collision 
process, starting from the collision between constituents and developing the parton shower, that 
increases the number of final-state particles by means of strongly ordered subsequent emissions. 
Eventually, the interface with a phenomenological hadronization model, enables the comparison 
with experimental data. For these reasons, they are routinely used by experimentalists to 
simulate signal and backgrounds processes in physics searches. Nevertheless, the demand for 
better and better predictions from high energy experiments calls for improving the precision 
of existing SMC's, including NLO corrections. The MCQNLO PQ method has shown first how to 
reach NLO accuracy for inclusive quantities, implementing the hard subprocess at NLO and 
developing showers within the leading logarithmic approximation, avoiding double counting of 
radiation. In this way one achieves benefits of both approaches: exclusive final states generation 
of SMC's and accuracy of NLO calculations. 

The POWHEG method is a different prescription for interfacing NLO calculations with parton 
shower generators. It was first suggested in ref. [2J, and was described in great detail in ref. [3]. 
This method is independent from the Monte Carlo program used for subsequent showering and 
generates positive weighted events only. In these respects it improves the MC0NL0 approach. 
Until now, the POWHEG method has been successfully applied to several processes, both at 
lepton (U [5] and hadron colliders [HJ [3 E HU1 HU [121 [131 [H]. In these implementations, it 
has been interfaced to the HERWIG [HJ[TB], PYTHIA [T7] and HERWIG++ [TBI SMC programs. 

In the POWHEG method the hardest radiatiorQ is generated first, independently from the 
following ones. SchematicalljH, the hardest radiation is distributed according to 



da = B ($b) d$B 



(i) 



x By hardest we mean the radiation with the highest transverse momentum, either with respect to the beam 
for initial state radiation (ISR), either with respect to another parton for final state radiation (FSR). 

2 Here we avoid entering into the details concerning the radiation regions and the correct treatment of the 
associated flavour configurations. The interested reader can find further explanations in the original papers (2]| 
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where B ($b) is the Born contribution and 



V(<S>! 



(2) 



is the NLO differential cross section at fixed underlying Born kinematics and integrated over 
the radiation variables. The transverse momentum of the emitted parton, with respect to the 
beam or to another particle, depending on the region of singularity, is denoted by Ut (^r)- The 
lower cutoff p™ m is necessary in order to avoid the coupling constant to reach unphysical values. 
V ($b) and R(&r) are the virtual and the real corrections and in the expression within the 
square bracket in Eq. ^ a procedure that takes care of the cancellation of soft and collinear 
singularities is understood, e.g. Frixione-Kunszt-Signer (FKS) Qjj] or Catani-Seymour (CS) 
dipolc subtraction. Then, 



exp 



ad 



B($ fl ) 



(fc T (® r )-Pt) 



(3) 



is the POWHEG Sudakov, that is the probability of not having an emission harder that pp. Equa- 
tion ([T]) can be seen as an improvement on the original SMC hardest-emission formula, since 
the Born cross section is replaced with B($b) which is normalized to NLO by construction. 
At small transverse momenta the POWHEG Sudakov becomes equal to a standard SMC one. 
However, the NLO accuracy of Eq. ([1]) is maintained for inclusive quantities. Moreover, the 
high— pt radiation region is correctly described by the real contributions 



da re B{<S> B . 



(4) 



since A7? re 1 and B/B re 1 + O (a s ). After having generated the hardest radiation, one can 
interface with any available shower generator, in order to develop the rest of the shower. To 
avoid the double-counting, the SMC is required to be either px— ordered or to have the ability 
to veto emissions with a pt harder than the first onq§. 



2 The POWHEG-BOX 

In a real collision process several colored massless partons are present, either in the initial or the 
final state. One thus should repeat the procedure outlined in Sec. [T]for every possible singular 
region, associated with any massless colored leg becoming collinear to another one, or soft. In 
order to do this, the full real emission cross section is decomposed into a sum of terms, each 
of which has at most one collinear and one soft singularities. The radiation is then generated 
independently in each of this regions, but only the hardest radiation is retained and the event is 
generated according to the flavour and kinematics associated to it. Because of this complexity, 
an automatic tool, the POWHEG-BOX, has been built [21], in order to help the inclusion of new 
processes. On the other hand, the POWHEG-BOX may also be seen as a library, where previously 
implemented processes are available in a common framework. 

The user wishing to include a new NLO calculation must only know how to communicate the 
needed information to the POWHEG-BOX. This happens either defining the appropriate variables, 
either providing the necessary Fortran routines. The required inputs are: 

3 All modern SMC generators compliant with the Les Houches Interface for User Processes |20j implement 
this feature. 
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1. The number of legs in Born process, e.g. nlegborn= 5 for pp — >■ [Z — > e + e ) + j. 

2. The list of Born and Real processes flavours flst_born(k=l . . nlegborn, j=l . . flstjiborn) , 
f 1st jreal (k=l . .nlegreal , j=l . . flstjireal) , accordin g to PDG convention^. Fla- 
vor is defined incoming (outgoing) for incoming (outgoing) fcrmion lines, e.g. [5 , 2 , 23 , 6 , 3 , 0] 
for bu — > Ztsg. 

3. The Born phase space Born_phsp(xborn) , given xborn(l . . .ndims) random numbers in 
the unit ndims-dimensional hypercube. This routine should also set the Born phase space 
Jacobian kn_jacborn, the Born momenta kn_pborn, kn_cmpborn in lab. and CM frames 
and the Bjorken x's, i.e. kn_xbl ,kn_xb2 . 

4. The routines that performs the initialization of the couplings, init_couplings, and the 
setting of the factorization and renormalization scales, set Jacjren_scales (muf ,mur) . 

5. The Born squared amplitudes B = |-M| 2 , summed and averaged over color and helici- 
ties as well as the color-ordered Born squared amplitudes Bjk and the helicity correlated 
Born squared amplitudes Bk^v, where k runs over all external gluons. The correspond- 
ing user routine is defined by setborn(p,bf lav, born, born j k , bornmunu) , for momenta 
p (0 : 3 , 1 :nlegborn) and flavour string bflav(l : nlegborn) . 

6. The real emission squared amplitudes 1Z user routine setreal(p,rf lav,amp2), for mo- 
menta p(0 : 3 , 1 :nlegreal) and flavour string rf lav(l :nlegreal) . 

7. The finite part of the interference of Born and virtual amplitude contributions Vb = 
2Re{£> x V}, after factorizing out the common factor TV = ^{i-e) (fp) ' ^ c user rou t mc 
is called setvirtual (p(0 : 3, 1 : nlegborn) , vf lav(l : legborn) .virtual) . 

8. The Born color structures in the large 7Y C limit, set by the routine borncolour_lh and 
passed through the Les Houches interface [20] . 

We remark that items (p~|[7|) are the usual ingredients needed to perform a NLO calculation 
in any subtraction method. Item (|8|) is instead needed to provide a defined color structure to 
the SMC generator. Internally, the P0WHEG-B0X implements the FKS subtraction procedure 
in a general way. At the beginning, it automatically evaluates the combinatorics, identifying 
all the singular regions and the corresponding underlying Born contributions. It also performs 
the projection of real contributions over the singular region and computes the subtraction 
countcrterms, from soft and collinear approximations of real emissions. Then, it builds the 
ISR and FSR phase space, according to the FKS parametrization of the singular region and 
performs the integration. Eventually, one gets the NLO differential cross section. At this stage, 
one can also interface to some analysis routine to obtain NLO differential distributions as a 
byproduct. After the integration stage, it performs the calculation of upper bounds for an 
efficient generation of Sudakov-suppressed events and then the generation of hardest radiation, 
according to the POWHEG Sudakov. At this point, the generated event, which contains at most 
only one extra radiation, has to be passed to a standard SMC program, for developing the rest 
of the shower. This can be done either on-the-fly or storing the events on a Les Houches events 
file [23j . Standard analysis routines, at partonic and hadronic level, are provided for included 

4 Internally gluons are labelled instead of the PDG value of 21. At the moment of writing the event on the 
Les Houches common block, the PDG value is restored. 
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processes, as well as drivers for common SMC generators. Users can modify them or implement 
new ones. 



2.1 Recent developments 

Recently, thanks to this framework, the relatively complex process of Z + 1 jet production 
has been implemented |24| . This is a promising processes for jet calibration with the early 
LHC data. It is also an important source of missing energy signal as well as a background to 
many new physics searches. In experimental studies carried up until now, the NLO theoretical 
calculations were supplemented by correction factors for shower, hadronization and underlying 
event effects. However, these factors were evaluated by means of standard LO SMC programs. 
It is clear the advantage to have a SMC program which is NLO accurate, in order to ease and 
improve the comparisons with experimental results. We have tried a simple approach [24] to 
merge consistently Z and Z + 1 jet samples, in order to obtain a description as smooth and 
accurate as possible both in the low and high-transverse momentum regions. The results are 
showed in Figs. [1] and [2] 




Figure 1: The pr distribution of the Z bo- 
son in single Z production (blue dashed 
curve), in Z + 1 jet (black solid) and in 
the merged sample (red solid). 



Figure 2: The pr distribution of the next- 
to-hardest jet in single Z production (blue 
dashed curve), in Z+l jet (black solid) and 
in the merged sample (red solid). 



From the two figures, one can see how the merged sample models both the single Z Sudakov 
form factor, that plays an important role in resumming collincar/soft logarithms in the low- 
Pt region and the high-p-r behaviour of the next-to-hardest jet, which follows the Z + 1 jet 
distribution. In this last figure, jest are reconstructed according to the kx algorithm, imposing 
also an angular separation Rjj > 3 . 
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